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We present a cosmological model in 1 + m + p dimensions, where in m-dimensional space there are 
uniformly distributed p-branes wrapping over the extra p-dimensions. We find that during cosmolog- 
ical evolution m-dimensional space expands with the exact power-law corresponding to pressureless 
matter while the extra p-dimensions contract. Adding matter, we also obtain solutions having the 
same property. We show that this might explain in a natural way why the extra dimensions are 
small compared to the observed three spatial directions. 
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I. INTRODUCTION 

The inflationary paradigm is successful in explaining 
away the basic shortcomings of standard cosmology, like 
the monopole, horizon and flatness problems. One might 
along this line claim that we now have a cosmological 
scenario which remains plausible very close to big bang. 
However, the possibility of an initial singularity still re- 
mains to be dealt with. One hopes that this problem 
will be resolved once we understand the theory of gravity 
down to sizes comparable to Planck length. String the- 
ory is one of the leading candidates for that but for its 
consistency one has to introduce extra dimensions. Ob- 
servationally if these extra dimensions exist their sizes 
are much smaller than the size of our perceived universe. 
Thus it is of importance to seek for cosmological models 
where this difference can be accommodated in a natural 
way. 

It is plausible that the sizes of all dimensions started 
out the same, possibly close to Planck length. After var- 
ious cosmological eras the perceived universe grew to its 
size we observe today. Therefore the problem is to ex- 
plain how the extra dimensions remained comparatively 
small. 

In what follows we propose a toy model, which we be- 
lieve gives an answer to this problem in a natural way. 
The basic motivation is a flavor of the idea once exposed 
by Brandenberder and Vafa Q (see also i and 0): when 
a p-brane wrap over the extra p-dimensions it resists ex- 
pansion, much like a rubber band would if wrapped and 
glued over the surface of a balloon. There has been con- 
siderable activity in the literature on similar ideas, see for 
example the subject of "brane gas cosmology" 

Our model is formulated in l-t-m-l-p-dimensions, where 
m and p refer to the observed and compact dimensions 
respectively, in which there is a uniform (with respect 
to m) distribution of p-branes wrapping over the extra 
p-dimensions. In order to incorporate this assumption it 
is enough to take the observed space to be topologically 
non-compact. At this point it is apparent that this ap- 
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proach is different than the brane gas cosmology, where 
it is postulated that branes can wrap anywhere. 

In this paper we focus on time dependent solutions 
to Einstein equations and flnd exact expressions where 
possible. In the earlier work on brane (or string) gas cos- 
mology the interest was mainly on the thermodynamical 
aspects and on the incorporation of T-duality invariance 
of string theory to cosmology (see for instance 0, 
and 0). For that reason, the problem was studied in 
the framework of dilaton gravity and it was found that 
the wrapped branes can prevent cosmological expansion 
of the internal space. As we will see, the same can be 
achieved in Einstein gravity without invoking T-duality 
invariance. 

Recently, dynamical aspects of brane gas cosmology 
in M-theory has been studied in ,15j and ^16]. Our ap- 
proach is technically similar. Namely, we also obtain the 
energy momentum tensor by coupling the brane action to 
the gravity action and assume a uniform distribution of 
such branes. However, we concentrate on pure Einstein 
gravity rather than M-theory. 

The main phenomenon we observed is that the ex- 
istence of wrapped p-branes makes the m-dimensional 
space grow in exactly the same way as pressureless mat- 
ter would, while the p-dimensional compact space is con- 
tracted following a power-law depending on the numeri- 
cal value of p. From this perspective the model has pre- 
dictive power since it gives a definite proportion between 
the sizes of the observed and the compact directions. 

We have also observed that adding ordinary matter 
does not change the mentioned behavior appreciably ex- 
cept for the case with negative pressure. For instance, 
adding radiation one still finds that the observed space 
expands and the compact space contracts. On the other 
hand, in case of vacuum domination both observed and 
compact spaces will grow exponentially. More on our line 
of reasonings is presented in Section III. 

The organization of the manuscript is as follows: In 
Section II we show how to wrap the p-branes over extra 
dimensions, present the resulting Einstein equations and 
the aforementioned solution. In Section HI we add ordi- 
nary matter to the energy- momentum tensor. In Section 
IV we give an estimate for the current size of the in- 
ternal dimensions using the solutions. In section V, we 
compare the relative strengths of the expansion and con- 
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traction forced by p-branes. The last section is devoted to 
conclusions and possible future extensions of the model. 



II. COSMOLOGY OF WRAPPED p-BRANES 



For more than one p-brane one has to change 6{x — Xq) 
with the sum '^niS{x — xi) where is the number of 
coincident p-branes at the position x\. We now assume 
that there are uniformly distributed such p-branes in Mm 
and do the following replacement 



Consider a D-dimensional space-time which has the 
following metric 

dt^ + e^^dx'dx' + e^^dy^dy", (1) 



ds^ 



where i — l,..,m, a = and the metric functions 

A, B, C depend only on time t. Here, is chosen to pa- 
rameterize the observed directions which we label as Mm 
and parameterizes the extra space labeled by Mp. We 
also define = {t,x'^,y°-). We would like to determine 
the cosmological evolution in the presence of p-branes 
wrapping over Mp. The dynamics of p-branes are deter- 
mined by the Polyakov action 
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(2) 

where = (r, cr") are world-volume coordinates, Tp is 
the p-brane tension, is the space-time and 7Q./3 is 
the world- volume metrics, and ^^(^) is the map from 
world-volume to space-time. It is easy to show that the 
following p-brane configuration in JQ) is an extremum of 
the Polyakov action (0) 



X 



T, 



(3) 



a 



a 



where Xq are constants giving the position of the p-brane 
in Mm- We assume that ?/" and a'^ coordinates are topo- 
logically so that there is no surface term coming from 
the variation of the action Q . We would like to calculate 
the back reaction of this configuration on the geometry. 
For that we couple (0) to Z?-dimensional Einstein-Hilbert 
action 



1 



d^Xyf^R, 



so that 



S — Se + Sp. 



(4) 



(5) 



The energy momentum tensor for the p-brane can be cal- 
culated from lO which gives 



d^p+il^j»PdaX''dpX''S{X-X{£_)). (6) 
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In the orthonormal frame {e^^dt, e^dx^, e'^dy'^), I© cor- 
responding to the p-brane configuration ^ takes the 
form 



T-.-. 
13 

T-i 

ab 



0, 



(7) 



-mB 



S{x - Xo)Sab- 



y^Tl; S{x 



Xl) 



dx' n{x') S{x — x'), (8) 



where n(x) is the number of p-branes per unit volume at 
the position x. Homogeneity of Mm implies that n{x) is 
constant and this gives the following energy momentum 
tensor 



Tit = i^Tpe 
T~'. = 

Tab = -nTpe 



-mB 



(9) 



-mB 



Sab 



One can easily check that V^T'"' = 0. We also note that 
a scaling of x coordinates x Xx requires, by definition, 
a scaling of n by n ^ n/A. 

In solving Einstein equations, we first impose the gauge 
choice A — mB+pC which fixes i-reparameterization in- 
variance in the metric In this gauge, the field equa- 
tions that follow from the action ((S)) can be written as 



A" - A'2 
B" 

C" 



mB'^ + pC'^ = C" 



m + p — 
m — 2 



m + p — 1 



where ' denotes differentiation with respect to t. The last 
two equations in (|10|) imply that B and C are propor- 
tional to each other upto the linear terms in t. Ignoring 
these terms B and C can be solved up to an undeter- 
mined integration constant. The first equation in H10() 
then fixes this integration constant. Switching to the 
proper time coordinate {which we again denote by t), we 
finally obtain the following metric 



ds^ = -dt^ 



where 



(at)-' 



-dx'dx' + (at) ^^ri+TTdy^d?/'', (11) 



n n T, 



TO^(p + 1)^ 
2(m + p — 1) (m — mp + 4p) 



p 



(12) 



For the physically important case of m — 3, from Hll|l . 
the scale factors R3 and Rp of the observed and compact 
dimensions respectively are determined as follows 



i?3(<) 

Rp{t) 



(at)-2/3(p+i) 



(13) 
(14) 



It is evident that the power-law of the observed space 
is exactly the same as the one for pressureless matter in 
standard cosmology. This is somewhat expected since the 
observed space part of energy momentum tensor vanishes 
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as it does for pressureless matter. On the other hand, 
(|14|l shows that even in pure Einstein gravity wrapped 
branes can prevent expansion of the internal dimensions. 
This is contrary to the general expectation in the litera- 
ture (see for instance , 0| and ) which is based on 
the intuition that negative pressure would increase the 
expansion rate (note that branes apply negative pressure 
along the wrapping directions) as in vacuum domination 
during inflation. However, H14(l indicates that negative 
pressure does not necessarily imply expansion in Ein- 
stein gravity. Moreover, even in de Sitter phase of the 
early universe, negative pressure would give exponential 
contraction with a negative Hubble constant. 

III. ADDING MATTER 

We now add ordinary matter to analyze a more realis- 
tic model. We will use for generality the following energy 



momentum tensor for matter 

Tjj^i, = diag(p, pi, Pa) , (15) 

where the indices refer to the obvious orthonormal frame 
in m and 

Pi = ^P, (16) 
Pa = vp, (17) 

where uj and v are constants. The energy-momentum 
conservation V^T^'' — determines p in terms of the 
metric functions 

p - Poe-('+"^'"''-('+"^'''', (18) 

where po is a constant and we again impose the gauge 
A = m,B +pC. In the presence of matter, Einstein equa- 
tions itTIH) are modified as follows 



A" - A'^ + mB'^ + pC'^ = C"-KVo(l + t^)e('""^™''+^'~"^^^, 

\p + l)nTp + po(l + (j5 - l)uj - pv)e-^'-B-{i+'^)pC 



{m+p-l)B" = K2gmB+2pC 
{m+p-l)C" = K2gmB+2pC 



(19) 



-(m - 2)nTp + po(l + {m - l)v - mw)e""'^-^-(i+'')'"^ 



These equations have a much richer solution space and 
it is not possible in this manuscript to exhaust every in- 
teresting one. There are various possibilities here de- 
pending on what one chooses for uj and u. However it is of 
crucial importance to accommodate for the inflationary 
paradigm, so we consider this first. In the 4-dimensional 
cosmology one chooses w = — 1 to achieve an exponential 
growth. This however also has a physical interpretation. 
The inflationary solution is the one in which vacuum en- 
ergy is dominant. It would be, we feel, quite unnatural 
to assume that a constant cosmological constant to be 
absent from the extra dimensions, so we belive it is nec- 
essary to chose V ^ —1 also. In this case, comparing 
the terms on the right hand side of (|19|l . it is clear that 
the functions multiplying po is times larger than the 
terms coming from p-brane sources. Evidently the for- 
mer term will dominate the equations in time granted 
the scale is increasing In this regime, p-brane 

sources can be self-consistently ignored (i.e. one can set 
Tp = 0) and the only contribution to the energy mo- 
mentum tensor comes from a cosmological constant, po, 
which will yield the usual exponential growth of the in- 
flationary period ^3 where all dimensions expand with 
the same exponent. 

To support the above argument, we also made a sim- 
ple numerical integration of the equations (|19|l in the 
proper time coordinate. We took p = 1, m = 3 and 



nTp/ Po — 10.0. Also the initial conditions for the run 
were such that near t = the scale factors obey (|13ll and 
(|14|l . The resulting plot is presented in Fig. 1. Note 
that the inflationary regime sets in around C" — 0. Be- 
fore this, the internal dimensions contract. Consequently, 
even though the e-foldings after the cosmological con- 
stant dominates are the same, the internal dimensions 
exit the inflationary period with a smaller final value for 
the scale factor than that of the observed dimensions. 
This difference depends on nTp/ p^. For our numerical 
run, the ratios of the scale factors (i.e. e^/e*-^) can be 
read from the graph and one finds that the scale factor of 
the internal dimensions is about 6.8 times smaller than 
that of the observed dimensions. 

In solving (|19|l for other cases of a; and v, we start with 
the following ansatz 

B = 6iln(t) +ln62, (20) 
C = ciln(i) +lnc2. (21) 

Using the last two equations in (|19|l one can uniquely 
determine these four constants and the first equation is 
satisfied identically. The constants which determine the 
power-law for expansion turn out to be 

2{\ + v) _ 2uj 

' " ~m{l-2LU + i^)' " p{l -2u: + v) ^^^^ 
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FIG. 1: Evolution of the scale factors B and C, for ui = v — 
— 1, p — 1, m — 3 and nTp/po — 10.0. See text for the initial 
conditions. 



Form (pn|) . and if^ . one can also fix the metric 

function A using the gauge A = mB + pC. Thus, all 
unknown functions in are now determined. Making a 
coordinate transformation to switch to the proper time 
coordinate {which we again denote by t), the metric can 
be written as 

ds^ = -dt^ + {ait)^dx'dx' + {a2t)'^^^dy°'dy'', (23) 

where ai and 012 depend on 62 and C2, and they can be 
set to 1 by scalings of x and y coordinates. 

It is clear from (I23f) that one should choose v ^ —1. 
Also, for the above ansatz to work, 62 and C2 should 
be determined to be positive numbers. This imposes, 
after a straightforward but somewhat tedious algebra, 
the following conditions on the parameters of the model 



{\ + v){\-v) uj{\-uj) 



> {I - LU + i^){lu - ly) {2A) 



(l + t^)(m-2) ^ ioip+l) 
m p 
l-2uj + v ^ 0. 



(25) 
(26) 



For pressureless dust uj ^ v ^ Q, and these conditions 
are satisfied for all m > 2 and p. In this case the scale 
factor for the extra dimensions is equal to 1; the expan- 
sion forced by the dust is compensated by the contraction 
forced by p-branes. 

For radiation, two different equations of state are pos- 
sible. When the sizes of all dimensions are close to each 
other one has oj ^ v ~ 1/ {m + p). In this case, (|24|) and 
(|26|l are satisfied identically. However, H25I) restricts the 
possible values of m and p. For instance, for m — Z l|25|) 
imposes that p > 1. On the other hand, when the extra 
dimensions are very small compared to the observed di- 
mensions one has oj = 1/m and v — Q. In this case (|24|) 
and (|26|l are satisfied identically, but (|25|) implies that 
mp > 3p -|- 1 and thus one should take m > 4. 



Before closing this section let us emphasize that these 
constraints only show that the ansatz chosen in solving 
the differential equations is not valid for all parameters. 
However, one can still infer a general dynamical behavior 
from the above solutions. We belive that for the cases 
where ()23|l is not valid, one would observe similar effects 
for ordinary matter coupled to p-branes. 



IV. AN ESTIMATE FOR THE SIZE OF THE 
INTERNAL DIMENSIONS 

Here we try to obtain an estimate for the current size 
of the extra dimensions along the lines of the model we 
presented. In this section we set to = 3. We take the 
universe to be filled with ordinary matter (character- 
ized with an equation of state depending on w, v) and 
wrapped branes. We further assume that after big bang 
all dimensions started out close to Planck length. The 
standard model of cosmology tells us that the universe 
passed through three different eras. First, an infiationary 
period took place. After inflation, there was a radiation 
dominated era followed by a matter dominated one which 
still (possibly T^) is going on. We assume that adding 
wrapped p-branes along the extra dimensions do not al- 
ter this history (but modify the power-law expansion as 
we discussed). 

During inflation, one should take lo — u — —1 which 
simply represents a positive cosmological constant. As we 
discussed in the previous section H19(l gives the usual ex- 
ponential growth (with a quantitative modification, see 
Fig. 1). Assuming approximately 70 c- foldings during 
inflation (which is required by cosmological phenomenol- 
ogy), the sizes of all dimensions grew to about 10^"''m 
from the Planck length. Note that we are not attempt- 
ing to answer questions like what is deriving inflation or 
how the graceful exit occurs. 

Just after the inflation during radiation dominated era 
one can set oj = v = l/(p-|-3). Note that the temperature 
at the beginning of the radiation era is expected to be of 
the order of IQ-'^^K which corresponds to a length scale 
10~^^m. Since this is much smaller than 10~^m, the radi- 
ation is allowed to apply pressure along the extra dimen- 
sions. In this period, the solution (|23|l can be used to de- 
scribe the cosmological evolution of the universe |23| . As- 
suming that the vacuum energy driving the inflation was 
about 10^"' GeV, the Hubble parameter of the inflation 



can be determined to be H^^ w 10""^^ sec. We match 
the solution H23I) to de Sitter space by demanding that 
the Hubble parameter of the observed universe is contin- 
uous. Therefore, H23|l is valid starting from t « 10^"^^ 
sec. Using that the radiation domination ended at about 
10^^ sec, ()23f) gives the size of the internal space at the 
end of radiation to be 10"^ x 10-i38/[p(p+4)]in. 

As the observed universe expands, pressureless matter 
started to dominate the cosmological evolution and one 
can then set lo = v — Q. During this period, (|23|l implies 
that the extra dimensions stay fixed. 
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Summarizing, we found that after an exponential 
growth during inflation, the extra dimensions contracted 
till the end of radiation dominated era and then re- 
mained unaltered. This gives the following estimate for 
the present size (vp) of the extra dimensions 

_ _ 138 

rp « 10"^ X 10 m. (27) 

For j3 = 1 we get ri « 10"'^^ m. For higher values of p one 
finds much larger estimates. However let us remind the 
reader of the digression we had on the structure of the 
inflationary scenarios. In reality, during the inflationary 
period the scale factor of the internal dimensions will 
be subject to a lesser growth than that of the observed 
dimensions. For instance, for the numerical case we have 
considered, rp should be scaled down by a factor of 6.8. 
Recall that this is for nTp/pQ — 10.0. Larger values for 
that ratio will make r„ smaller. 



V. MORE ON p-BRANE COSMOLOGY 

In section we found that in the presence of uni- 
formly distributed wrapped p-branes, the observed space 
expands while the compact space contracts. In this sec- 
tion, we would like to compare relative strengths of these 



two effects by assuming existence of two branes where 
they will be uniformly distributed wherever they don't 
wrap. Let us consider the following metric in 1+m+p+q- 
dimcnsions 

+ e^^dz^^dz''^ (28) 

where ai,5i — l,...,p, 02,62 — l,---,? and the metric 
functions depend only on time t. We assume that in 
there are p and q dimensional branes wrapping over y 
and z coordinates, respectively. Following our reason- 
ings in section^ one can easily write down the energy 
momentum tensor corresponding to this configuration 

rp rp —mB — qD 1 rp — mB—pC 

-'-it ~ -'p^ -'q^ 5 

^ 0, 

T,,i„ = -T,e-'"^-^^, (29) 

T-i ^ rp —niB—pC 

0262 ~ "J ' 

where Tp and Tq are respective brane tensions and hatted 
indices refer to the orthonormal basis in (|28|l . At this 
point one can check that V pT^'^ = 0. Imposing the gauge 
A = mB -j-pC -j- qD the Einstein equations can be written 
as 



A" - A'^ + mB'^ + pC'^ + qD'^ = C" + D" - B", 

{m+p + q- l)B" = «;2gmi3+pC+gD ^ l)npTpeP^ -t- (q + l)n,T,e«^] , 

{m+p + q- 1)C" = «;2gmi3+pC+«D ^ ^ _ 2)npTpeP^ + (q + l)n,r,e«^] , (30) 

{m+p + q -1)D" = K^e"'''+P^+^''[ip+l)npTpeP^ -im+p-2)nqTqe'^''], 



where Up and Uq are number of branes per co-moving 
volumes parameterized by (a:, z) and {x, y) coordinates 
respectively. To solve these equations exactly we first 
note that the last three equations in (|30|l imply that 



(m - 3)B + ip+ l)C +iq+l)D = 



(31) 



up to the terms linear in t which we ignore in the follow- 
ing. Using 131|) . the middle two equations can be solely 
expressed in terms of B and C. To find a power-law 
solution we choose 



B = bi \n{t) +ln62, 
C^ci ln(t)+lnc2. 



(32) 
(33) 



It is now straightforward to check that the four constants 
can be determined uniquely. For instance we find 



bi 



Cl 



-2(p- 



2pq) 



9pq + ■m{p + q — pq) ' 

2(m- 3)g 
'dpq + m{p + q - pq) 



(34) 
(35) 



Using H31|) and the gauge choice A = mB+pC + qD, one 
can also determine metric functions A and D. The first 
equation in H30|l is then to be checked for consistency. 
A straightforward calculation shows that it is satisfied 
identically. Let us also emphasize that, contrary to the 
case we encountered in section III, the constants 62 and 
C2 turn out to be positive for all possible values of to, p 
and q. 

Using the above results, in the proper time coordinate 
the metric can be written as 

ds^ = -dt^+Rl^dx'dx'+Rldy'''dy'''+Rldz''^dy''' (36) 

where 

1 ^ 2(p + g + 2pq) 
ln(i?„) = ■ — ■ ■ In(ami) 



ln(i?p) 
h-i{Rq) 



2>pq + m{p + q + pq) 
2q{m~3) 
3pq + m{p + q + pq) 

2p{m — 3) 
ipq + m{p + q + pq) 



\n{apt) (37) 
ln{aqt) 
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and am, ctp and Uq are positive constants. Comparing 
(|S7|l with the results of section II, we see that wrapping 
a p-brane and a g-brane is physically distinct from wrap- 
ping a p + g-brane. The main difference is that in the 
former case the p-branes are uniformly distributed along 
the dimensions over which g-branes wrapped and vice 
versa. 

Another interesting aspect is that only the dimension 
of the observed space m determines whether the compact 
dimensions expand or contract. For the two brane con- 
figuration considered above we see that m = 3 is the crit- 
ical value where the internal dimensions are stabilized by 
the cancellation of the expansion forced by the uniformly 
distributed p-branes with the contraction forced by the 
wrapping of g-branes. We find that the compact dimen- 
sions expand for m > 3 and they diminish for m < 3. 

Generalizing and (|37|l . one would guess that the 
exponent of the power-law of the internal dimensions is 
proportional to m — fc — 1 where k is the number of par- 
titionings of the extra dimensions, i.e number of distinct 
brane configurations. For large fc, one should increase m 
to have contracting extra dimensions. 

VI. CONCLUSIONS AND FUTURE 
DIRECTIONS 

We have shown that if one allows for p-branes wrapping 
around extra dimensions Einstein equations allow for so- 
lutions where their size diminishes during cosmological 
evolution. In this paper, we have focused on the main as- 



pects of the idea and omitted various details which may 
form topics for future studies. 

First, since the p-branes allow for power-law solutions 
for the observed universe which are exactly the same as 
pressureless matter one can think of them as a form of 
dark matter. 

Second, it is an important question to ask about the 
time of forming of wrapped p-branes during cosmological 
evolution. In this work we have assumed that they did 
exist since the beginning of time but it is plausible that 
they may form later. One important question along this 
line is about the connection of the wrapped branes to the 
exit from inflation. 

Let us also recall that Brandenberger and Vafa argued 
that string interactions may yield an upper bound 
(to < 3) on the number of observed dimensions if one 
assumes that all the directions are compact. We note 
in this context that the two brane scenario we have dis- 
cussed in section V yields a lower bound for to if one 
insists on preventing expansion of the extra dimensions 
during entire cosmological history. One is tempted to 
speculate that the two ideas might be merged to fix to. 

Finally, one could further think of adding ordinary 
matter to the configuration discussed in section V. This is 
expected to make the internal dimensions grow for to > 3. 
Of course a detailed numerical analysis might prove oth- 
erwise. One further extension could be to wrap branes 
over branes. In this case we expect the effects to add 
up in the directions the branes intersect. And this might 
result in a hierarchy of scale factors. 
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Note that both the cosmological constant and wrapped 
p-branes force B to increase. This argument would fail 
if, for instance, we would have the function C instead of 
B since the cosmological constant and the p-branes have 
opposite effects on C. 

To see this more explicitely, set Tp = in 11911 . Then, 
it is easy to show that A = — ln{Ht) and B — C = 
— ln{Ht)/{m -|- p) is a solution where — po{m + 
p)/(m-|-p— 1). Making a further coordinate transforma- 
tion, one can see that (0 represents de Sitter space. 
Recent observations indicate that the universe is accel- 
erating which suggests that the energy density is now 
dominated with some kind of undetermined energy (dark 
energy) having negative pressure. In this paper we are not 
going to discuss possible modifications implied by the ex- 
istence of dark energy to our scenario. 
Since in 12311 the internal space contracts, the equation 
of state for radiation also starts to shift. As the extra 
dimensions become much smaller than the scale set by 
the temperature, one should set o) = 1/3 and u = Q. We 
ignore this subtlety in the following since we try to obtain 
a rough estimate to see if the scenario makes sense in a 
first crude approximation. 



